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In a digital fiber optical communication system, the optical power re- 
quired at the receiver input to achieve a desired error rate depends upon 
the shape of the received pulses. In systems employing multimode fibers 
and/or broadband sources, we can experience pulse spreading in propaga- 
tion because of the group velocity differences of different modes or because 
of dispersion. In an effort to control or compensate for pulse spreading, 
we can trade off coupling efficiency between the light source and the fiber 
(by varying the core-cladding index difference or bandlimiting the source), 
scattering loss in the fiber (by introducing mode coupling), and equalization 
in the receiver at baseband. This paper investigates the optimal trade-off 
for various fiber-source combinations. 

I. INTRODUCTION AND REVIEW OF BACKGROUND MATERIAL 

In digital fiber optic communication systems, as in other digital 
systems, the received power required at a repeater to achieve a de- 
sired error rate depends upon the shape of the received pulses. A 
previous paper 1 showed that the minimum average power requirement 
results from a pulse that is sufficiently narrow so that its energy spec- 
trum is almost constant for all frequencies passed by the receiver 
(ideally, an impulse). For other received pulse shapes, we can define 
the additional power required, in decibels, as a "power penalty" for 
not having impulse-shaped pulses. Typical calculations of this power 
penalty for "on-off" signaling and a receiver employing avalanche 
gain with a high impedance front end 1 are shown for various families 
of received pulse shapes in Fig. 1. In that figure, the parameter <t/T 
is defined as follows : 

i = t> { z l h *w m - [ j /mo**] 1 } , (i) 
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Fig. 1 — Typical power penalty vs <r/T. 

!T = spacing in time between binary digits 
h p (t) = received optical pulse shape 
A = area under h p (t). 

We shall refer to a- as the rms pulse width. * 

It has been shown 2 - 3 that, in long fibers, the "power impulse re- 
sponse" of the fiber approaches a Gaussian shape. In the rest of this 
paper, we assume that the received optical pulse is Gaussian in shape 
and that it has an rms width determined by the fiber delay distortion. 
That is, we assume that the rms width of the fiber input pulse is 
sufficiently small so that the power penalty associated with the rms 



* In a Gaussian-shaped pulse, the rms width is about 0.425 times the full width 
between half-amplitude points. In a rectangular pulse, the rms width is 1/V12 the 
full width. 
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sum of that width and the rms fiber impulse response width is the same 
as the power penalty associated with the rms fiber impulse response 
width alone. 

From various heuristic analyses (see the appendix), we can conclude 
that the rms width of the received pulse is approximately the rms sum 
of the delay distortion in the fiber resulting from material dispersion 
(because of the variation of group velocity with wavelength associated 
with the use of a broadband source) and the delay distortion associated 
with the spread in the group delays of various fiber modes (when a 
multimode fiber is used). That is, 

O" = (^dispersion T ""mode) y (2) 

where «T dl8pei . 8lon ~ optical bandwidth • fiber length = BL and where 
©"mode is determined as follows 4 

Case 1. Conventional clad multimode fibers without mode coupling. 

A/7 

°mode = 0.289 L, 

c 
where 

n = index of refraction of the core 

A = (index of refraction of the core — index of refraction of the 

cladding) /n 
c = speed of light. 

Case 2. Conventional clad fibers with complete mode coupling after 
a distance L c - 

(r mode = 0.289 — JUTc for L > L c , 
c 

= 0.289 — L for L ^ L c . 

c 

Case 3. Ideal graded index fiber without mode coupling. 

°mode = 0.037 L. 

c 

Case 4. Other fibers 2 can be treated once the techniques outlined below 
are understood. 

From the definitions of o- dl8perslon and <r mode above, we see that, for a 
broadband (incoherent) source, the material dispersion contribution 
to o- can be controlled by limiting the optical bandwidth B being used. 
However, if the optical source bandwidth, B , must be reduced by 
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filtering, then the average power into the guide will be reduced by the 
factor B/B . Similarly, we can control the mode delay spread by reduc- 
ing the index difference A. If a multimode (incoherent) source is being 
used, the average power into the guide is proportional to A. Thus, we 
trade off input power against mode delay spread. Furthermore, even 
when we use a coherent source which in principle can be focused into 
any fiber, we must be careful when A becomes significantly smaller than 
0.005, since the fiber loss at bends becomes large. (Exactly what value 
of A is too small to be practical is an open question.) For fibers with 
mode coupling, we can control <r mo de by decreasing Lc (increasing the 
mode coupling). However, this causes the coupling to radiating modes 
to increase, thus increasing the fiber loss. 5 Here again, there is a trade- 
off between the average power we receive at the fiber output and the 
received rms pulse width. For a fixed shape of the mechanical spectrum 
of fiber geometry perturbations, the radiation loss per unit length of 
the fiber resulting from mode coupling is inversely proportional to Lc, 
i.e., 

radiation loss in nepers = a L/Lc, (3) 

where 

a = constant depending upon the shape of the mechanical spectrum 
of the geometry perturbations causing coupling (and possibly 
upon the index difference A). Lc depends upon the amplitude of 
the mechanical perturbations. 

In the following sections we derive the optimal trade-off between a, 
B, A, and Lc for various combinations of sources and fibers to maximize 
the allowable fiber length L between the optical source and the 
repeater. 

II. ANALYSIS 

2.1 Incoherent source, conventional clad fiber, no mode coupling 

Let the average power into the guide be P a when the index difference 
A is at some maximum practical value A and when the full source 
optical bandwidth B is being used. Let the loss of the fiber be a nepers 
per kilometer. Let the power penalty from the nonzero value of the 
received rms pulse width, in nepers, be j{a/T). Let the required power 
at the receiver be P r when <t/T = 0. If we use a value of A ^ A and 
filter the source output to have an optical bandwidth B ^ B , then we 
must have 

P '-TW e ~" L ^ ?* MT \ (4) 

A ti 
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where, from (2), 

'=|( c, k l ) !+ ( c ^ l )T =w+ ^' 

and C\, C 2 are constants. 

We rewrite (4) as follows (using equality to maximize L), 

P T 6 ~ \A B e 

To maximize L, we must choose A/A and B/B to maximize the term 
in braces subject to the constraint that these ratios cannot exceed 
unity. We define —10 log (term in braces) as the "excess loss." 

By equating appropriate partial derivatives of the excess loss to 
zero, we obtain the following equations for optimizing A and B (A and B 
which minimize excess loss). 



(5a) 



''(t)^ 1 


provided A/A„ ^ 1, 




otherwise A/A = 1, 


1'(-t)£^ 


provided £/B < 1, 




otherwise B/B„ = 1, 



(5b) 



where f'(z) = d/dx^f(x)^\ x=z and <r m and a d are defined in (4). For 
sufficiently long lengths L, where both A/A and B/B are less than 
unity, we obtain [by adding (5a) to 5(b)] 

dm = Td = ~J2 > W 

where x is the solution f'(x)x/2 = 1. More specifically, we obtain the 

following 

r BL _ xT AL _xT m 

Cl B~ -^' ° 2 A7~ v2 {n 

and therefore* 

[x^T 2 1 

for 

L _ , ( x x | 

^ maximum of j ^ , ^ j ■ 

* Throughout this paper we use the parameter L/T (guide length/time slot width) 
frequently. The larger the fiber length or the smaller the time slot width, the more 
excess loss must be incurred to control or compensate for pulse spreading. 
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From the Gaussian power penalty curve of Fig. 1, we obtain the value 
of x where f'(x)x = 2 to be 0.37. At that value of x, 

- 10 log e-' (x) = 3.3 dB. 

As L/T decreases, either A/A or J3/B will eventually reach unity. 
When that happens, a will approach a* or a m , respectively, for shorter 
lengths L. Then, from (5), a/T will approach the solution of f'(z)z = 1. 
Furthermore, the excess loss will approach either 



excess loss 



10 log 



[&*"] 



(8) 



(9) 



if A/A reaches unity first and L/T > z/Ci or 

excess loss — > — 10 log jfT e~ f{z) 

if 6/ B reaches unity first and L/T > z/C 2 . 

From the Gaussian curve of Fig. 1, the solution of f(z) z = 1 is 
z = 0.3. At that value of z, -10 log e~ fM = 1.8 dB. 
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Fig. 2 — Excess loss vs L/T for conventional fiber. 
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For L/T sufficiently small so that both A/A and B/B = 1, the 
excess loss will approach zero as L/T decreases. 

Example 1 

In a fused silica fiber, the material dispersion has been measured at 
9 ps/km per angstrom of wavelength difference. With a typical GaAs 
LED* (light-emitting diode), this results in a value of ad of 1.5 ns/km 
of pulse spreading at the full bandwidth B . Thus, C\ can be set at 1.5 
ns/km. For a fiber with a maximum index difference A of 0.01, C% 
would be given by 14.5 ns/km. 

From (7) we obtain 

excess loss = -10 log 2 J^%^ m5) <r'» W 
= -20 log ^ + 28.3 dB 

for L/T > 0.174, where the length L is in kilometers and the time slot 
width T is in nanoseconds. 

From (8) we find that, for 0.0206 < L/T < 0.174, the excess loss 
asymptotically approaches 



excess loss — > — 10 log 



»1 = -10 log ^ + 18.6 dB 



LC 2 



for 0.0206 < L/T < 0.174. 

For L/T < 0.0206, the excess loss asymptotically approaches zero. 
Figure 2 is a plot of excess loss vs L/T in this example. 

2.2 Incoherent source, ideal graded index fiber, no mode coupling 

Following the same procedures as hi 2.1, we can replace o m for a 
self-focusing fiber by C 3 (A/A )\L (where A„ is the maximum allowable 
value of A). We then obtain the following set of equations which 
determine the values of A and B that minimize the excess loss 

''(i)^- 1 prided i/A„gi, (9a) 

otherwise A/A = 1, 
= 1 provided £/B c < 1, ,„, -, 

otherwise B/B„ = 1- 



' \T)aT 



* Assuming a Gaussian-shaped optical spectrum with bandwidth between the half- 
power points of about 400 A. 

DIGITAL FIBER OPTIC COMMUNICATION 791 



For sufficiently long lengths L, where both A/A„ and 8/B„ are less than 
unity, we obtain [by adding (9a) to twice (9b)] 



^1 
v3 






(10) 



<?d — Cijj- L — x 



excess loss = — 10 log 



where x' is the solution of f'(x')x' = 1.5, and where we must have 



L . , I x 

7= > maximum ot < -~ 



^vf and kyjl] 



For the Gaussian power penalty curve shown in Fig. 1, we have 
x' = 0.34 and -10 log e~^ x,) = 2.6 dB. 

As before, as L/T decreases, either £/B or A/A will reach unity. 
Thereafter, for smaller values of L/T we have the following: either 



zT 

a — > <r d ; excess loss — > — 10 log j^- e~ /i ' ) , 



(11) 
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Fig. 3 — Excess loss vs L/T for graded index fiber. 
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where f'(z)z = 1, provided A/A reaches unity first and L/T > z/Ci, or 

o- — > <r m ; excess loss — ► — 10 log I I =77- j e-/<*'> , 

where f'(z')2z' = 1, provided £/£ reaches unity first and L/T > z'/C 3 . 
For the Gaussian power penalty curve of Fig. 1, we obtain z = 0.3, 
-101oge-' ( *> = 1.8 dB, z' = 0.24, -10 log g-M = 0.98 dB. 

When L/T is sufficiently small so that A/A„ and B/B are both equal 
to unity, then the excess loss asymptotically approaches zero for 
smaller L/T. 

Example 2 

From Example 1 we have Ci typically 1.5 ns/km. From (2) we have 
C 3 typically 0.019 ns/km for an ideal graded index fiber with a maxi- 
mum A = A of 0.01. 

From (10) we obtain 

excess loss = — 15 log ( y J + 4.85 dB 

for L/T > 10.3, where L is in kilometers and T is in nanoseconds. 
From (11) we obtain 

excess loss -*■ — 10 log ( 7 ) + 8.78 dB 



(1) 



for 0.2 < L/T < 10.3. 

For L/T < 0.2, the excess loss asymptotically approaches zero as 
L/T approaches zero. 

Figure 3 shows a plot of excess loss vs L/T for this example. 

2.3 Incoherent source, conventional fiber with adjustable mode coupling 

Now consider a conventional fiber with adjustable mode coupling. 
Using the same notation as in 2.1, we have the following condition 
from (2), (3), and (4)' 

A R 

p e -aL e -a LIL c ^ P^fKolT)^ ( 12 ) 



*As mentioned in Section I, the parameter a depends upon the shape of the 
mechanical coupling spectrum and possibly upon the index difference A. Since this 
dependence upon A is not known analytically, we assume a = constant independent 
of A. One could also assume a <* A* for some iV (probably negative) and still obtain 
simple results sinuTar to those that follow using analogous techniques. 
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where 



-((c.|*) , + (<?.Evn5) , ) , -w+*p 1 



provided L > Lc. 

By setting appropriate partial derivatives to zero, we can minimize 
the excess loss given by 



excess loss = — 10 log 



A B 



The optimizing values of A, fi, and L c satisfy the following equations 
(we are assuming a fixed by the shape of mechanical coupling 
spectrum) . 

''(I)!- 1 ' P^vided i/A„g 1, (Ua) 

otherwise A/A = 1, 

>'(?);l = 1 ' provided B/B ° - *• a4b) 

otherwise $/.B = 1, 

xf = /' ( f ) ^ > P rovided L > Zc ' (14c) 

^ = 0.5, provided A/A ^ 1 and L > L c - (14d) 

For sufficiently long fibers and if a ^ 0.5, we will have L > L c , 
A/A < 1, £/Bo < 1, and therefore the following will hold:* 

ex eeas loss =-10 log [ 2 -gg^e— ], (15) 

where x is the solution of f'(x)x/2 = 1 

provided 



«oL _ n K a sr £ ^r_ 

X;~ U - 5 ' A 2C 2 LV^' 5 vXLC, 



It is convenient to consider L/T and L/Z c as separate parameters. 

* It is interesting to note that, with optimal mode couplin g, in th e region where 
A < A , the optimal value of A is increased by the factor V0.5/a„ relative to the 
no-mode coupling case [see formulas for A/A„ in (7) and (15) and also (8) and (16)]. 
Further in this region (A < A ), the excess radiation loss from mode coupling (a L/L e ) 
is always 0.5 neper. 
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As L/T decreases, either A/A or 11/ B will reach unity. 
If B/B reaches unity first, then the excess loss will asymptotically 
approach the following for smaller values of L/T. 



excess loss — ► — 10 log 



e -/u) e -o.5i (16) 



C 2 L-y/2^ 



-r- = 0-5\ 



lc ' ' A C 2 LV2^ ' 

where f'(z)z = 1, provided 6/B reaches unity first and 

L/T > z/(C 2 V2^,). 

If A/A reaches unity first, then the excess loss will asymptotically 
approach the following for smaller values of L/T 



excess loss — * — 10 log 

2? ^ zT 
B LCt ' 



zT ,. , 
LCx 



(17) 



provided A/A reaches unity first and L/T ^ z/C\. 

For values L/T below that at which A/A and B/B both equal unity, 
the excess loss asymptotically approaches zero. 

Example 3 

Using the same parameter values as in Example 1 and assuming* 

a„ = 0.1, we obtain the following 

T 
excess loss = -20 log j + 27 dB 

for L/T > 0.174 km/ns, 

T 
excess loss => — 10 log y + 17.3 dB 

for 0.0462 < L/T < 0.174, and 

excess loss => for ™ < 0.0462. 

Figure 4 shows a plot of excess loss vs L/T for this example. 



* At this point, the achievable value of o in practical fibers is a subject of specula- 
tion. We choose a = 0.1 arbitrarily. 
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Fig. 4 — Excess loss vs L/T for conventional fiber with coupling (a = 0.1). 

2.4 Laser source, conventional fiber with adjustable mode coupling 

Here we assume that, to avoid excessive loss at bends, the index 
difference in the fiber, A, is fixed at some minimum allowable value 
Amin. To control pulse spreading we can trade off mode-coupling radia- 
tion loss against equalization penalty. The condition we must satisfy is 

p e -aL e -a LILc ^- p e /ClT) 

We wish to choose L c to minimize the excess loss given by 

excess loss = -10 log {«— 'M-ce-'W), (18) 

where* 

a = Ci^lLLc, C 4 = 0.289 A^n/c. (19) 

We obtain the optimizing equation : 

«.t/i c = g/'(^)f- (20) 



* Material dispersion is assumed negligible for a coherent source. That is, we assume 
a single mode and a short-term bandwidth less than 1 A. 
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Fig. 5 — Excess loss vs L/T — laser source, conventional fiber A m i„ = 0.001, o = 0.1. 

To solve (20) we can pick a value of cr/T and solve for f'{a/T) and 
f((r/T) graphically from Fig. 1. We then use those results in (20) 
to solve for L/L c . Then we substitute into (19) to find L/T and into 
(18) to find the total excess loss. 

Example 4 

Using Amm = 0.001 and a = 0.1, Fig. 5 shows a plot of excess loss 
vs L/T for this example. 

111. APPLICATIONS 

If the optical power required at the receiver when the received pulses 
are very narrow is P r and the transmitted power (at maximum band- 
width and index difference) is P, and if the fiber loss in the absence of 
mode coupling loss is aL, then we must have 

10 log (P.e- aL ) - excess loss (L) = 10 log (P r ) 

or, equivalently, ' 

p t 

10 log y e ~ aL = "excess gain (L)" ^ excess loss( L). 



* We define excess gain as number of decibels by which P.e~ aL exceeds P r . In effect, 
it is equal to the "allowable excess loss." 
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At a given bit rate, l/T, and given a, P., and P r , we can plot excess 
loss (L) and "excess gain (L)" simultaneously. The intersection of the 
two curves gives the maximum allowable distance L between the 
transmitter and the receiver. 

Example 5 

Assume that, at a bit rate of 25 Mb/s {T = 40 ns), the required 
received power P r is approximately —58 dBm. Assume that a conven- 
tional fiber with mode coupling (a„ = 0.1) and loss a = 5 dB/km is 
used. Assume that an incoherent source is being used with P„ = —13 
dBm for A max = 0.01. Assume that Ci and C 2 are 1.5 and 14.5 ns/km 
so that Fig. 4 applies. Figure 6 shows a plot of excess loss and excess 
gain vs L. It is apparent that the maximum length L between the trans- 
mitter and the receiver is 6.8 km. At that distance, the excess loss 
= 11.5 dB. Further, we have A S 0.0024, £/B = 1 and L c = 1.36 
km.* 

IV. COMMENTS AND CONCLUSIONS 

The purpose of this paper has been to show how we can combine 
analytical results on fibers and repeaters to determine maximum re- 
peater spacings by optimization of available parameters. Since the 
fiber art is still young, many assumptions above are subject to ques- 
tion. We can summarize a few possible criticisms here. 

It is not known whether the assumption of the Gaussian pulse shape 
leads to overly conservative estimates of the equalization penalty. 
With time and experiments, the Gaussian pulse shape approximation 
will probably be improved upon. 

It is not known yet how much control the designer will have over 
the mode coupling and the index difference. Future analyses will 
have to take into account the practical constraints on these parameters. 

It is not known yet whether optical filters of the type assumed above 
can be built. Further, the above analysis neglects in-band insertion 
loss. 

It is hoped that, although the above analysis is somewhat simplistic, 
it can serve as a guide to the fiber system designer by pointing out the 
concepts and trade-offs involved. 



• It is interesting to note that, if mode coupling were not allowed, the excess loss 
curve would be about the same (calculated from Fig. 2) and therefore the maximum 
length L would be the same. However, at the maximum length A would be 0.0011. 

798 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1 974 



10 



20 



25 



30 



EXCESS GAIN 



I I I I 




J I I I I I I 



0.1 



100 



L IN KM 

Fig. 6 — Excess loss and gain vs L. 
APPENDIX 

We wish to show heuristically that the total rms width of the fiber 
impulse response is the rms sum of the contribution from dispersion 
and the contribution resulting from mode delay spread. 

Suppose that, if the fiber is excited by a narrow-band source at wave- 
length X, the resultant output response is h\(t). Let the mean arrival 
time and rms width of h\ (t) be defined as 



T\ 



<x\ 



= -^ fthx(()dt 

'{if 



(21) 



Pfa(t)dt - tI 



where 



A\ = I h\(t)dt = area of h\(t). 



Now suppose the fiber is excited by a narrow pulse from a broadband 
source having its output distributed in wavelength according to the 
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spectrum S(\). Intuitively,* we can write the fiber output response 
as follows: 

h(t) = js(\)h x (t)d\ = jlS(\)A,2 (^f) dk. (22) 

Define 5(X) as S(X)A X and let 

5 = fs(\)Aid\ = fh(t)dt. 

Let 5(X) = lS{\)Ax/ST\. It follows from (22) that the rms width of 
h{t) is given by 

-(j/*»-(s/*®)T 

= |[/^(\) f /x] 

+ I" JT*(\)S{\)d\ - (jT(X)S(X)dxYj| 4 . (23) 

The first term in square brackets in (23) is a weighted average of the 
mean square width of the narrow-band pulse at different wavelengths. 
The second term is the mean square deviation of the narrow-band- 
mean-arrival time, i.e., the dispersion a\. If we next assume that 
a\ « constant = a m (i.e., that the rms width of the narrow-band im- 
pulse response is not dependent upon wavelength within the band of 
interest), then we obtain 

<T= {oi + 02}*, 

which is the desired result. 
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* But not rigorously. This is where the argument becomes heuristic. 
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